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We report low-temperature thermal conductivity k of pyrochlore Yb 2 Ti 207 , which contains frus¬ 
trated spin-ice correlations with significant quantum fluctuations. In the disordered spin-liquid 
regime, k(H) exhibits a nonmonotonic magnetic field dependence, which is well explained by the 
strong spin-phonon scattering and quantum monopole excitations. We show that the excitation en¬ 
ergy of quantum monopoles is strongly suppressed from that of dispersionless classical monopoles. 

Moreover, in stark contrast to the diffusive classical monopoles, the quantum monopoles have a very 
long mean free path. We infer that the quantum monopole is a novel heavy particle, presumably 
boson, which is highly mobile in a three-dimensional spin liquid. 


Rare-earth pyrochlore oxides exhibit various exotic 
magnetic properties owing to their strong geometrical 
frustration experienced by coupled magnetic moments 
on the tetrahedral lattice (Fig. 1(a)) [1]. The most ex¬ 
plored materials are Ho 2 Ti 2 07 and Dy 2 Ti 2 07 , in which 
the magnetic moments can be regarded as classical spins 
with a strong easy-axis (Ising) anisotropy [1,2]. The frus¬ 
tration of these moments results in a remarkable spin ice 
with macroscopically degenerate ground states, in which 
each tetrahedron has the “two spins in, two spins out (2- 
in-2-out)” configuration. This spin structure is character¬ 
ized by dipolar spin correlations with a power-law decay, 
which is observable as the unusual pinch-point shape of 
spin structure factor by neutron scattering [3, 4]. One of 
the most remarkable features of the spin-ice state is that 
it hosts emergent magnetic monopole excitations; the 
first excitation is 3-in-1-out configuration [5, 6]. This pro¬ 
duces a bound pair of north and south poles, which can 
be fractionalized into two free magnetic monopoles. This 
classical monopole excitations are gapped and dispersion¬ 
less (Fig. 1(b)). Therefore the propagation of monopoles 
occurs only diffusively and the monopole population de¬ 
cays exponentially at temperatures well below the gap. 
Of particular interest is how the spin-ice ground state 
is altered by the quantum fluctuations, which may lift 
the degeneracy of the spin-ice manifold, leading to a new 
ground state such as quantum spin-ice state [7-11]. To 
clarify this issue, uncovering newly emergent elementary 
excitations in the presence of quantum fluctuations is 
crucially important. Although exotic excitations such as 
gapless photon-like mode have been proposed theoreti¬ 
cally, the nature of excitations are poorly explored. 

Among the magnetic pyrochlore materials, Yb 2 Ti 2 07 , 
Er 2 Ti 2 07 , Pr 2 Sn 2 07 and possibly Tb 2 Ti 2 07 host strong 
transverse quantum fluctuations of magnetic dipoles 


owing to speudospin-1/2 of magnetic rare earth ele¬ 
ments [12-15]. In particular, Yb 2 Ti 2 07 with well sep¬ 
arated crystalline electric-field excited levels is a good 
model system to study the influence of the quantum 
effects on monopole excitations [16]. Yb 2 Ti 2 C >7 under¬ 
goes a weakly first-order ferromagnetic phase transition 
at Tq ~ 0.2 K [17-20]. It is widely believed that quan¬ 
tum fluctuations melt the spin-ice state into a spin-liquid 
state. Therefore, the pinch point structure observed by 
neutron scattering above Tq indicates the presence of a 
spin-liquid phase with spin-ice correlations [17]. The full 
set of Hamiltonian parameters is determined by inelastic 
neutron scattering experiments, providing a prototypical 
system described by an effective pseudospin-1/2 quantum 
spin-ice model [12]. The Hamiltonian consists of three 
main interactions, Jn, J± and J z ±. Here J\\(= 2K) is 
the Ising component of the nearest neighbor interaction, 
J_l(= 0.58K) is the XY-component and J z ±{= 1.7K) is 
the off-diagonal component. Finite J± and J z ± produce 
quantum fluctuations (Fig. 1(c)). 

Here, to study the elementary excitations in the spin 
liquid state of Yb 2 Ti 2 07 , we measured the thermal con¬ 
ductivity, which is a powerful probe for low energy exci¬ 
tations at low temperatures, providing a sensitive mea¬ 
surement of a flow of entropy conducted by magnetic ex¬ 
citations and phonons. The thermal conductivity has 
been reported in the classical spin-ice state of Dy 2 Ti 2 C >7 
recently. However, the interpretation of the thermal con¬ 
ductivity of Dy 2 Ti 2 07 appears to be complicated owing 
to the strongly suppressed phonon thermal conductiv¬ 
ity by unknown additional scatterings (see Supplemental 
Material [21]). In fact, suggested heat transport by clas¬ 
sical monopole is at odds with the diffusive motion of 
the dispersionless classical monopoles. We show that the 
thermal conductivity of Yb 2 Ti 2 07 is rather simple: the 
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FIG. 1. (color online), (a) Spin-ice structure on frustrated py- 
rochlore lattice, (b) Magnetic monopole excitations in classi¬ 
cal and quantum spin ice. In classical spin ice the gap energy 
is twice the Ising interaction of magnetic moments 2 Jn. In the 
quantum spin ice the off-diagonal interaction J z ± gives rise 
to a dispersive monopole excitation. The photon excitations 
based on the XT-component J± lift the ice degeneracy in 
the ground state, (c) Collective motion of quantum magnetic 
monopoles. 

phonon term shows a B/T scaling and the monopole con¬ 
tribution vanishes below Tq as expected. Our analysis 
shows the evidence of the substantial heat transport by 
quantum magnetic monoples, whose excitation energy is 
significantly suppressed from that of classical monopoles. 
The quantum magnetic monopoles are highly mobile due 
to quantum fluctuations, in stark contrast to the local¬ 
ized and diffusive nature of classical monopoles. 

High quality single crystals of Yb 2 Ti 207 were grown 
by the floating zone method. Thermal conductivity was 
measured along [1 ,-1,0] direction by the standard steady- 
state method in a dilution refrigerator. Magnetic field 
was applied along [1,1,1] and [0,0,1], perpendicular to the 
heat current. Specific heat was determined by the quasi- 
adiabatic heat pulse method in a dilution refrigerator. 

Figure 2(a) shows the temperature dependence of ther¬ 
mal conductivity divided by temperature k/T in zero 
field and at /iqH = 12 T measured on a single crystal 
of Yb 2 Ti 207 . Distinct jump in k/T at zero field is ob¬ 
served at Tq. As shown in Fig. 2(b), the specific heat 
C of the single crystal taken from the same batch shows 
a sharp and large jump at Tq [17]. We note that in 
the previous studies, such a sharp single jump in C/T 
had been reported only in the powered samples [19, 22], 
demonstrating the high quality of the present crystal. 
We also note that the pinch point features in neutron 
scattering has been clearly observed in the single crystal 
which shows a similar specific heat jump [17]. As shown 
in Fig. 2(c), zero-field k/T above Tq follows a T-linear 
dependence with negligibly small intercept at T = OK. 
The absence of residual k/T |t^ok in the spin-liquid 
state with spin-ice correlations will be discussed later. 

As clearly seen in Fig. 2(a), magnetic field strongly en¬ 
hances the thermal conductivity. Figures 3(a) and 3(b) 
and their insets show the field dependence of k{H)/T for 



FIG. 2. (color online), (a) k/T at zero and /iqH — 12T 
applied along [1,1,1] direction is plotted against temperature. 
The heat current is applied along [1,-1,0]. At Tc, k(0)/T 
exhibits a jump, indicated by an arrow. Inset illustrates the 
measurement configuration of the thermal conductivity, (b) 
Specific heat divided by temperature C/T at zero field, (c) 
k/T at zero field in the spin liquid state above 0.2 K. Grey 
line is a fit to a T-linear dependence k/T = AT with A = 
0.15 W/K 2 m. 


different field directions. As illustrated in the inset of 
Fig. 3(c), there are three characteristic regimes; low-field 
regime (i) where k(H)/T decreases with i7, intermediate 
field regime (ii) where k(H)/T increases, and high-field 
regime where k(H)/T exhibits a saturation. 

In the present system, heat is transferred by phonons 
and magnetic excitations: k = k p + K m . We point out 
that the field dependence of k(H)/T in the (ii) and (iii) 
regimes are dominated by the phonon contribution k v de¬ 
termined by spin-phonon scattering, which contains elas¬ 
tic and inelastic processes. The elastic scattering (de¬ 
termined by J\\) is enhanced with increasing disorder of 
spins and thus this scattering process should be mono- 
tonically suppressed by the alignment of spins with in¬ 
creasing magnetic field. A recent calculations of mag¬ 
netoresistance in a fluctuating spin-ice state indicates 
that the electron-spin elastic scattering rate decreases 
with increasing magnetization [25], which supports this 
trend. The inelastic scattering is directly related to the 
quantum dynamics of spin. In this inelastic scattering, 
the leading spin-flip process accompanies a hopping of 
a monopole to the neighboring tetrahedron (which is re¬ 
lated to Jj_), because this process requires much lower en¬ 
ergies than creation or annihilation of monopoles. This 
scattering is suppressed with field by the formation of 
Zeeman gap. Therefore an external magnetic field sup¬ 
presses both elastic and inelastic scatterings, leading to 
the enhancement of the phonon thermal conductivity k p . 

In the regime (iii), the Zeeman splitting energy g^sH 
well exceeds both of the magnetic interactions and ther¬ 
mal energy, ggsH Jy, Jl, J z ± and &bT. In this situa- 
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tion, where all spins are fully polarized and the magnetic 
(spin-wave) excitations are gapped with a gap gusH, 
thermal conductivity is almost entirely dominated by 
the pure phonon contribution without spin scattering be¬ 
cause of the following reasons. First, elastic spin-phonon 
scattering is absent due to the perfect alignment of spins. 
Second, inelastic scattering is also absent due to the for¬ 
mation of the large Zeeman gap. Third, spins do not 
carry the heat due to the Zeeman gap. Since purely 
phononic thermal conductivity is insensitive to magnetic 
field, k(H)/T in the regime (iii) is nearly independent 
of H. In the regime (ii), the phonon mean free path is 
significantly reduced by the spin-phonon scattering due 
to the spins thermally excited across the Zeeman gap. 
In fact, as shown in Fig. 3(c) which plots k/T as a func¬ 
tion of gsH/ksT , all data collapse into a single curve 
except for the low tiuH/ksT regime. The fact that data 
for both field directions stabilizing different spin config¬ 
urations (3-in-l-out for H || [1,1,1] and 2-in-2-out for 
H || [0, 0,1]) follow the same curve implies that the elas¬ 
tic spin-phonon scattering dominates over the inelastic 
scattering in this regime. It is intriguing that the H/T 
scaling curve appears to follow the Brillouin function (the 
dashed line in Fig. 3(c)). Here we fitted the data with 
^-factor of 0.79, which is comparable to Lande ^-factor 
(8/7) of Yb. This coincidence with the Brillouin function 
calls for further theoretical investigations. 

A particularly important information for the elemen¬ 
tary excitations is provided by k(H)/T in the regime (i), 
where k(H)/T decreases with H (the insets of Figs. 3(a) 
and (b)) and exhibits striking deviations from the H/T- 
scaling curve (Fig. 3(c)). This low-field behavior of 
k(H)/T is most likely due to the monopole contribution 
because of the following reasons. First, the initial 
reduction with H cannot be explained by spin-phonon 
scattering, which always increases k(H) with H as dis¬ 
cussed above. Second, the deviations in regime (i) appear 
below T* ^ 4K, where the pinch point features in neu¬ 
tron scattering appears [17]. In addition, T* is close to 
the temperature 2 J||/&b, above which monopole excita¬ 
tion disappears. Third, as shown in the inset of Fig. 3(a) 
(see Supplemental Material [21] for more detail), the ini¬ 
tial reduction of k(H)/T disappears below Tc, which is 
consistent with the monopole scenario because ferromag¬ 
netic ordering prevents the monopole formation. 

The decrease of k(H)/T with H implies that the num¬ 
ber of monopoles is reduced with H at low fields. This 
reduction is expected even in the dispersionless classi¬ 
cal monopoles with gap 2Jy (see Supplemental Mate¬ 
rial [21]). However, in the classical case, the number 
of monopoles will decay exponentially with decreasing 
temperature below T* ^ 2J^/ks • Therefore the ob¬ 
served quite substantial reduction of k(H)/T even at low 
temperatures well below T* is totally inconsistent with 
the classical monopoles. The results provide strong evi¬ 
dence that the monopole excitation gap is dramatically 



FIG. 3. (color online), (a) Field dependence of k/T of 
Yb 2 Ti 2 07 for H || [1,1,1] with the heat current along [1,- 
1,0]. The inset shows k(H)/T at low field. Solid red and 
black lines are k/T(H) at 0.18 and 1 K, respectively. Data are 
shifted vertically for clarity, (b) The same plot for H || [0,0,1]. 
Double-headed red arrow in the inset of (a) indicates the ini¬ 
tial reduction of k(H)/T for H || [0,0,1] at T=0.6K, which is 
estimated to be 0.03 W/K 2 m, giving a lower-bound estimate 
of the monopole contribution, (c) Normalized thermal con¬ 
ductivity ft/ftsat plotted against /xbH/&bT, where ft S at is the 
saturated thermal conductivity at high fields. ft S at at high 
temperatures is determined so as to fit the scaling curve. The 
dashed line represents the Brillouin function with spin=l/2, 
assuming #=0.79. The inset illustrates the typical behavior of 
k(H)/T. There are three characteristic field regimes, (i),(ii) 
and (iii) indicated by different colors. 


suppressed from the classical monopole, suggesting the 
emergence of dispersive quantum magnetic monopoles il¬ 
lustrated in Figs. 1(c) and (d). We also note that the 
substantial reduction of monopole density by the low 
field will result in a reduction of the inelastic spin-phonon 
scattering process related to the monopole hopping dis¬ 
cussed above, which further emphasizes the significant 
role of the quantum monopoles themselves as a heat con¬ 
ducting carrier at low fields. 

As shown in the inset of Fig. 4, k(H) decreases as 
k(H) = /c(0) — olH 2 (a > 0) at very low fields. As 
the thermal conduction by magnetic excitations is de¬ 
termined by the number of low-energy itinerant quasi¬ 
particles, this is a measure of the suppression rate of 
magnetic monopoles at low fields. Figure 4 depicts the 
temperature dependence of a for H || [0,0,1] and [1,1,1]. 
As the temperature is lowered, a first increases, decreases 
after showing maximum at T max = 0.3-0.5K and sud¬ 
denly vanishes at Tc- The difference in the magnitude 
of a in the two field directions may be related to the 
expected difference in the density of 3-in-l-out configu- 
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FIG. 4. (color online). Temperature dependence of the initial 
slope of k(H) determined by fitting k(H) — a(0) — aH 2 , for 
H || [0,0,1] and [1,1,1]. Inset shows k(H) plotted as a function 
of H 2 at very low field. 


ration at high fields, but the trends of a(T) are similar 
in both cases. The enhancement of a with decreasing 
T below T* can be accounted for by the reduction of 
thermal smearing of the monopole excitations. The re¬ 
duction of a at lower temperatures is expected when the 
thermal energy scale ^T max becomes comparable to a 
fraction of the monopole excitation gap (see Supplemen¬ 
tal Material [21]). In addition, possible ferromagnetic 
fluctuations, because of the weakly first-order nature of 
the transition at Xc, would also suppress a near Tq. In 
any case, the observed low energy scale of fer max indi¬ 
cates that the excitation gap is strongly reduced than the 
classical monopole case (2Jy = 4K). The present results 
lead us to conclude that the thermally excited quantum 
monopoles carry substantial portion of the heat partic¬ 
ularly in the regime (i). This is reinforced by the fact 
that k/T at zero field shows a distinct decrease below Tq 
(Fig. 2), where the phonon contribution k v is expected to 
be enhanced owing to the ferromagnetic spin alignment. 

Next we demonstrate that quantum monopoles are 
highly mobile in the crystal lattice. Assuming the ki¬ 
netic approximation, the monopole contribution to the 
thermal conductivity is written as K m = C m v£/ 3, 
where C m is the monopole contribution in the specific 
heat, v is the velocity and £ is the mean free path of the 
monopoles. We estimate £ at 0.6 K simply by assuming 
that the amount of initial reduction of k(H)/T shown red 
double-headed arrow in the inset of Fig. 3(b) is attributed 
to the monopole contribution. The total specific heat 
C ~ 1 J/Yb-molK at 0.6 K (the inset of Fig. 2) and v, 
which is roughly determined by v ~ aJ z ±/2nh ~ 15 m/s, 
where a(= 0.43 nm) is the distance between neighbor¬ 
ing tetrahedra, yield £ ~ 100 nm, or equivalently the 
scattering time r ~ 2.5 ns. We stress that this long £ 


is still underestimated, since the total specific heat and 
the initial reduction of thermal conductivity give only 
an overestimate and underestimate, respectively, for the 
monopole contribution. This indicates that the excita¬ 
tions are mobile to a very long distance, £ > 250a, with¬ 
out being scattered. We stress that £ is much longer 
than the inter-monopole distance, which is estimated to 
be at most 5a, assuming monopole density of 1% of total 
number of tetrahedra. This corresponds to a very large 
coherent volume including more than ^ 10 7 tetrahedra, 
demonstrating highly mobile transport of this long-lived 
particle, whose effective mass is as heavy as ^ 2000 times 
the bare electron mass [26]. This extremely small scat¬ 
tering rate may be due to the quantum feature, which 
prohibits the simple monopole-antimonopole pair anni¬ 
hilation that violates energy conservation. 

The present results indicate that the quantum fluctu¬ 
ations dramatically change the nature of the monopole 
excitations. We note that the observed nearly gapless 
excitations are not relevant to the “photon” excitations 
predicted by Ref. [7-11]. This is because the character¬ 
istic photon energy E vhoton ~ ~ 0.05 K is one 

order of magnitude smaller than the present tempera¬ 
ture range, and hence the strongly temperature depen¬ 
dent a is incompatible with the photon excitations. The 
highly mobile heavy quantum monopoles in the spin liq¬ 
uid state is the most salient feature of the elementary ex¬ 
citations in frustrated magnetic pyrochlore systems with 
strong quantum fluctuations. Nearly ballistic propaga¬ 
tion phenomena of fractionalized magnetic excitations 
in spin-liquid states have been reported in spin-1/2 ID 
Heisenberg chain [27, 28] and 2D triangular lattice with 
antiferromagnetic interactions [29]. In the former ele¬ 
mentary excitation is spinon which obeys semion statis¬ 
tics [30] and in the latter excitation has been discussed in 
terms of spinon which obeys fermionic statistics [31-36]. 
In the present 3D system elementary excitation in the 
spin liquid state is quantum monopole, which is another 
fractionalized spinon. The residual k/T |t-^ok, which 
is distinctly present in the 2D case [29, 36], is absent 
in Yb 2 Ti 207 (Fig. 2(c)), implying that this 3D spinon 
is unlikely to be fermionic. In fact, bosonic spinon has 
been presumed theoretically in 3D pyrochlore lattice [37]. 
In ID Heisenberg system, the mean free path is infinite 
at nonzero temperature due to the integrability of the 
Hamiltonian. The highly mobile fermionic spinons in 2D 
and bosonic quantum monopoles in 3D may be a key 
feature of the elementary excitations in highly frustrated 
quantum magnets and its origin is an open question. 

We thank L. Balents, K. Behnia, S. Fujimoto, 
H. Kawamura, S. Onoda, and K.Totsuka for useful dis¬ 
cussions. Financial support for this work was provided 
by Grants-in-Aid for Scientific Research from the Japan 
Society for the Promotion of Science (JSPS). 












5 


[1] S. T. Bramwell, and M. J. P. Gingras, Science 294 , 1495 

( 2001 ). 

[2] A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan, 
and B. S. Shastry, Nature 399 , 333 (1999). 

[3] S. T. Bramwell, M. J. Harris, B. C. den Hertog, M. J. P. 
Gingras, J. S. Gardner, D. F. McMorrow, A. R. Wildes, 
A. L. Cornelius, J. D. M. Champion, R. G. Melko, and T. 
Fennell, Phys. Rev. Lett. 87 , 047205 (2001). 

[4] T. Fennell, P. P. Deen, A. R. Wildes, K. Schmalzl, D. 
Prabhakaran, A. T. Boothroyd, R. J. Aldus, D. F. Mc¬ 
Morrow, and S. T. Bramwell, Science 326 , 415 (2009). 

[5] C. Castelnovo, R. Moessner, and S. Sondhi, Nature 451, 
42 (2008). 

[6] D. J. P. Morris, D. A. Tennant, S. A. Grigera, B. Klemke, 
C. Castelnovo, R. Moessner, C. Czternasty, M. Meissner, 
K. C. Rule, J.-U. Hoffmann, K. Kiefer, S. Gerischer, D. 
Slobinsky and R.S. Perry, Science 326 , 411 (2009). 

[7] N. Shannon, O. Sikora, F. Pollmann, K. Penc, and P. 
Fulde, Phys. Rev. Lett. 108 , 067204 (2012). 

[8] O. Benton, O. Sikora, and N. Shannon, Phys. Rev. B 86, 
075154 (2012). 

[9] M. J. P. Gingras, and P. A. McClarty, Rep. Prog. Phys. 
77 , 056501 (2014). 

[10] L. Savary, and L. Balents, Phys. Rev. Lett. 108 , 037202 

( 2012 ). 

[11] M. Hermele, and M. P. A. Fisher, L. Balents, Phys. Rev. 
B 69 , 064404 (2004). 

[12] K. A. Ross, L. Savary, B. D. Gaulin, and L. Balents, 
Phys. Rev. X 1, 021002 (2011). 

[13] R. Applegate, N. R. Hayre, and R. R. P. Singh, T. Lin, 
A. G. R. Day, and M. J. P. Gingras, Phys. Rev. Lett. 109 , 
097205 (2012). 

[14] M. Hirschberger, J. W. Krizan, R. J. Cava, and N. P. 
Ong, arxiv: 1502.02006 (2015). 

[15] J. S. Gardner, M. J. P. Gingras, and J. E. Greedan, Rev. 
Mod. Phys. 82 , 53 (2010). 

[16] J. A. Hodges, P Bonville, A Forget, M Rams, K Krolas, 
and G Dhalenne, J. Phys. Cond. Matt. 13, 9301 (2001). 

[17] L.-J. Chang S. Onoda, Y. Su, Y. -J. Kao, K. -D. Tsuei, 
Y. Yasui, K. Kakurai and M. R. Lees, Nat. Commun. 3, 
992 (2012). 

[18] Y. Yasui, M.Soda, S. Iikubo, M. Ito, M. Sato, N. Ham- 
aguch, T. Matsushita, N. Wada, T. Takeuchi, N. Aso, and 
K. Kakurai, J. Phys. Soc. Jpn. 72 , 3014 (2003). 

[19] K. A. Ross, L. R. Yaraskavitch, M. Laver, J. S. Gardner, 
J. A. Quilliam, S. Meng, J. B. Kycia, D. K. Singh, Th. 
Proffen, H. A. Dabkowska, and B. D. Gaulin, Phys. Rev. 
B 84, 174442 (2011). 

[20] J. A. Hodges, P. Bonville, A. Forget, A. Yaouanc, P. 
Dalmas de Reotier, G. Andre, M. Rams, K. Krolas, C. 
Ritter, P. C. M. Gubbens, C. T. Kaiser, P. J. C. King, 
and C. Baines, Phys. Rev. Lett. 88, 077204 (2002). 

[21] See Supplemental Material for additional phonon scatter¬ 
ing in Dy 2 Ti 2 07 , disappearance of initial decrease of k(H) 
below Tc and theoretical field dependence of monopole 
density for classical spin ice. 

[22] A. Yaouanc, P. Dalmas de Reotier, C. Marin, and V. 
Glazkov, Phys. Rev. B 84, 172408 (2011). 

[23] G. Kolland, O. Breunig, M. Valldor,M. Hiertz, J. Friel- 
ingsdorf, and T. Lorenz, Phys. Rev. B 86, 060402 (2012). 

[24] G. Kolland , M. Valldor, M. Hiertz, J. Frielingsdorf, and 


T. Lorenz, Phys. Rev. B 88, 054406 (2013). 

[25] M. Udagawa (unpublished) 

[26] L. Pan, N. J. Laurita, Kate A. Ross, Edwin Kermarrec, 
Bruce D. Gaulin, and N. P. Armitage, arXiv: 1501.05638 
(2015). 

[27] A. V. Sologubenko, K. Gianno, H. R. Ott, U. Ammerahl, 
and A. Revcolevschi, Phys. Rev. Lett. 84, 2714 (2000). 

[28] K. Kudo. S Ishikawa, T. Noji, T. Adachi, Y. Koike, K. 
Maki, S. Tsuji, and K. Kumagai, J. Phys. Soc. Jpn. 70, 
437 (2001). 

[29] M. Yamashita, N. Nakata, Y. Senshu, M. Nagata, H. 
M. Yamamoto, R. Kato, T. Shibauchi, and Y. Matsuda, 
Science 328, 1246 (2010). 

[30] F. D. M. Haldane, Phys. Rev. Lett. 67, 937 (1991). 

[31] L. Balents, Nature 464, 199-208 (2010). 

[32] S.-S. Lee, and P. A. Lee, Phys. Rev. Lett. 95, 036403 
(2005). 

[33] S.-S. Lee, P. A. Lee, and T. Senthil, Phys. Rev. Lett. 98, 
067006 (2007). 

[34] M. S. Block, D. N. Sheng, O. I. Motrunich, and M. P. A. 
Fisher, Phys. Rev. Lett. 106, 157202 (2011). 

[35] M. Barkeshli, H. Yao, and S. A. Kivelson, Phys. Rev. B 
87, 140402 (2013). 

[36] D. Watanabe, M. Yamashita, S. Tonegawa, Y. Oshima, 
H.M. Yamamoto, R. Kato, I. Sheikin, K. Behnia, T. 
Terashima, S. Uji, T. Shibauchi and Y. Matsuda, Nat. 
Commun. 3, 1090 (2012). 

[37] T. Senthil, Discussion Meeting: Quan¬ 

tum entanglement in macroscopic matter, 
(htt p: / / www. ict s. res. in/lecture / details /1642/). 



1 


Supplemental Materials for ’’Thermal 
conductivity of quantum magnetic 
monopoles in the frustrated pyrochlore 
Yb 2 Ti 2 0 7 ” 


ADDITIONAL PHONON SCATTERINGS IN 

Dy2Ti207 

k/T of Y 2 Ti 2 0 7 , Dy 2 Ti 2 0 7 at zero field and Yb 2 Ti 2 0 7 
at H =12 T are compared in Fig.SI. k/T of Dy 2 Ti 2 0 7 at 
high temperature is strongly suppressed from Y 2 Ti 2 0 7 . 
It has been reported that k/T of Dy 2 Ti 2 0 7 decreases 
monotonically with H , indicating even smaller phonon 
k/T [SI]. This implies the presence of unknown phonon 
scatterings, which are not likely caused by crystal field 
excitations because of the large energy gap of first excited 
state, ^380 K [S2]. As discussed in the main text, k/T 
of Yb 2 Ti 2 0 7 at H =12 T represents the purely phononic 
k/T at low temperatures, where k^T « /jl-qH. k/T at 
2 K is much larger than Dy 2 Ti 2 0 7 and close to Y 2 Ti 2 0 7 , 
suggesting the absence of such unknown scatterings. 



FIG. SI. (color online). k/T of Y 2 Ti 207 , Dy 2 Ti 207 at 
zero field and Yb 2 Ti 2 C >7 at H —12 T applied parallel to [1,1,1] 
with the same heat current direction along [1,-1,0]. Data of 
Y 2 Ti 207 , Dy 2 Ti 207 are taken from Ref. [SI]. 


DISAPPEARANCE OF INITIAL REDUCTION 
OF k(H)/T WITH H BELOW T c =0.19K 

In a temperature range of 0.2 K< T <0.3 K, a clear 
kink appears in k(H)/T for H || [1,1,1] at the transition 
between spin liquid and ferromagnetic states, indicated 
by arrows in Fig. S2. The kink is shifted to lower field 
with decreasing temperature and vanishes below Tq in 
the ferromagnetic state. The position of the kink agrees 
with the field dependence of Tq(H) determined by spe¬ 
cific heat measurements, which increases with H at low 



FIG. S2. (color online). Field dependence of k/T of 
Yb 2 Ti 207 for H || [1,1,1] with the heat current along [1,-1,0]. 
Arrows indicate the field-induced transition from paramag¬ 
netic to ferromagnetic states. 



FIG. S3. (color online). Field dependence of k/T of 
Yb 2 Ti 207 for H || [0,0,1] with the heat current along [1,-1,0] 

field region [S3]. 

As discussed in the main text, the initial reduction 
of k(H)/T with H observed above Tq indicates thermal 
conduction of magnetic quantum monopoles. The initial 
reduction in the temperature range 0.2 K< T <0.3 K is 
interrupted by the field-induced ferromagnetic ordering 
and k(H)/T shows a characteristic enhancement with 
the field in the ferromagnetic state. The enhancement 
is understood by the suppression of elastic scattering of 
phonon due to the ordering of magnetic moments. As 
the kink disappears in the ferromagnetic state, the ini¬ 
tial reduction, which is the signature of monopole heat 
conduction, also disappears, in consistent with the sup- 
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TABLE SI. (color online). The energy of spin configurations 
in a single tetrahedron 

pression of spin-ice correlations below Tq reported by 
neutron scattering experiments [S4]. The initial reduc¬ 
tion disappears below Tq also for H || [0,0,1](Fig. S3). 

INITIAL REDUCTION OF CLASSICAL 
MONOPOLE DENSITY WITH MAGNETIC 
FIELD 

By calculating the classical monopole density (/? 3 i, 3- 
in-1-out and l-in-3-out configurations) in magnetic field, 
we show that p%i decreases with H 2 at zero-field limit, 
regardless of the field direction. 

Hamiltonian of a nearest-neighbor spin ice model is 
written as 

H = J Y, 

<i,j> j 

where J, which corresponds to Jy in the main text, is 
the nearest neighbor Ising interaction, the spin Sj is an 
Ising spin: Sj = crjdj, crj = ± 1 , with the anisotropy 
axes, d 0 = [1,1, l]/>/3, di = [1,-1,-1 ]/a/3, d 2 = 
[-1,1,-1]/V3, d 3 = [—1, —1, l]/\/3- With h = gn B H, 
the energy of all the spin configurations are shown in 
Table SI. 

For H || [1,1,1], p 3 i is derived as, 



FIG. S4. (color online). Classical monopole density p%\ 
is plotted against normalized magnetic field gpsH/J for 
H ||[1,1,1] and [0,0,1]. Dotted vertical line indicates the field 
gpBH/J= 3, at which level crossing of 3-in-l-out and 2-in-2- 
out occurs when H is applied along [1,1,1]. 



FIG. S5. (color online). h 2 coefficient of p 3 i, 031 , is plotted 
against normalized temperature, ksT / J. 
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8 cosh 3 ( 2hj/3t ) 
^[ 1 , 1 , 1 ] 


^[ 0 , 0 , 1 ] = 2 exp (2/£) [2 + cosh (4/ij/ V3t)] 

+ 8 cosh ( 2 hj/y/?>t) + 2 exp (— 6 / 1 ) 


^[ 1 , 1 , 1 ] = 6exp(2/t)cosh(4/ij/3t) 

+ 8cosh 3 (2foj/3£) + 2 exp (— 6 /t) 

Here, hj = h/J and t = ksT/J are normalized field and 
temperature, respectively. 

For H ||[0,0,1], 

8 cosh ( 2hj/\fZt ) 


The resulting field dependencies at different tempera¬ 
tures are plotted in Fig.S4. For H || [0,0,1], the energy of 
one 2-in-2-out configuration decreases the most by Zee- 
man effect, leading to monotonic increase of 2 -in- 2 -out 
density. As a result, P 31 decreases monotonically. For 
H || [1,1,1], the energy of one 3-in-l-out configuration de¬ 
creases the most and crosses with the lowest energy of 
2-in-2-out configuration at gpBH/ J =3. As this crossing 
occurs, P 31 rapidly increases with H. On the other hand, 
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in low field region, p 3 i decreases with field isotropically. 
We verify this by expanding p 3 1 with hj around hj= 0. 
For both the two field directions, hj-linear term vanishes 
and h 2 j term is identical. The isotropic field dependence 
of psi at zero-field limit is then, 

Pziihj) = P3i(0) - «31 hj 

_ 8 (exp ( 2 /t) — exp (— 6 /t)) 

31 3£ 2 (exp (— 6 /t) + 3 exp (2/t) + 4 ) 2 

The h 2 j coefficient a^\ is plotted against t = ksT/J 
in Fig.S5. It exhibits a maximum at kBT maXia31 / J =0.8. 
The relation between T max ^ a31 and the monopole excita¬ 
tion energy A 31 is then Asi=2J=2.bkBT rnax ^ a31 . In the 
main text, the initial H 2 decrease of k/T is ascribed to 
decreasing number of monopoles and the H 2 coefficient 
a in the field dependence of k/T exhibits a maximum 
at T max =0. 3-0.5 K. If T max is related to the monopole 


gap energy, it corresponds to 0.75-1.25K of monopole ex¬ 
citation energy, which is strongly suppressed from the 
classical one, 2Jy=4K [S5]. It should be noted, however, 
that this estimation is of purely classical monopole. 
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